Quantum engineering of photon states with entangled atomic ensembles 
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We propose and analyze a new method to produce single and entangled photons which does not 
require cavities. It relies on the collective enhancement of light emission as a consequence of the 
presence of entanglement in atomic ensembles. Light emission is triggered by a laser pulse, and 
therefore our scheme is deterministic. Furthermore, it allows one to produce a variety of photonic 
entangled states by first preparing certain atomic states using simple sequences of quantum gates. 
We analyze the feasibility of our scheme, and particularize it to: ions in linear traps, atoms in optical 
lattices, and in cells at room temperature. 

PACS numbers: PACS 



The deterministic generation of coUimated single and 
entangled photons is of crucial importance in Quantum 
Information, like in quantum cryptography [l|, quantum 
computation [2L quantum lithography Q or quantum 
interferometry J, • Most of the methods tested so far 
require high-Q cavities, something which is very demand- 
ing in practice [6, J,, Se M ■ The engineering of quantum 
states in atomic systems is now possible thanks to the 
experimental progress experienced by the field of Atomic 
Physics during the last years. In fact, with trapped ions 
it has been already possible to create so-called W [lOj 
and GHZ [T2| states of up to 8 ions. At the same time, 
scientists have been able to produce other kinds of entan- 
gled states with atoms in optical lattices. Further- 
more, with the advent of Rydberg techniques [13] it will 
soon be possible to create W-like states in that system or 
in atomic ensembles at room temperature. Apart from 
their fundamental interest, some of those states may have 
applications in precision spectroscopy [3, • 

In this work we show that the ability of creating those 
atomic states may have a strong impact in different sub- 
fields of quantum information, as it may lead to a very 
efficient way of creating certain kind of entangled pho- 
tonic states which are required in various applications. 
The main idea is to use a laser and an internal level con- 
figuration such that we can map the atomic state onto 
photonic states corresponding to modes propagating in a 
well defined direction. Our scheme uses the well known 
fact [T^ITtI that, under certain circumstances, light scat- 
tering takes place predominantly in the forward direction 
due to an interference effect. In fact, this effect is the ba- 
sis of one of the building blocks of the repeater scheme 
proposed in [l8| , and has been recently demonstrated in 
a series of experiments [l^ US [2l|. There, a single ex- 
citation is created in an atomic ensemble by detecting a 
photon emission in a certain direction. Then, the excita- 
tion is released in the forward direction by using a laser. 
Building on this fact, we propose to create certain kind 
of excitations by using quantum gates or atomic interac- 
tions, which give rise to the desired entangled states when 
they are released using a laser, and which propagate in 
the desired direction due to the mentioned interference 
effect. 



Let us consider a set of N atoms with (ground) hype- 
fine levels \g) and \sa,b) (see Fig. [T] (a)). We consider 
states of the form 

and linear combinations thereof. Here, |0) — \g)i. . . \g)N, 
and 
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where j excites an atom from \g)j to \sx)j, and are 
the equilibrium position of the atoms. In the limit 
N, Eq. Il]) defines a set of orthonormal collective states 
with Tlx atoms excited in \sx) and linear momentum k^. 
Those states can be indeed readily created using trapped 
ions or Rydberg techniques (see Appendix IAII3P . 

In order to release the photons, one sends a laser pulse 
of wavevector kL which couples level \sx) to some elec- 
tronically excited ones le^), respectively. The large pop- 
ulation of level \g) together with the initial entangle- 
ment (coherences) between the atoms, will now stimu- 
late the emission of photons from the excited states to 
the level \g), which overall will produce the mapping be- 
tween these states and the photonic states, 
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that is, ([T]) is mapped to a Fock state of photons with 
momenta k^ + kL and polarization a^, where is the 
polarization of the light in each decay channel. More- 
over, due to the linearity of this process, superpositions 
of states of the form ([T]) will be mapped onto superpo- 
sitions of photonic states Q. For example, the atomic 
state (|k(i),q(i)) + |q(i),k(i))) /V2 wih emit a pair of 
entangled photons in different directions. The mapping 
([3]) is strictly valid under ideal conditions, and in the limit 
N ^ oo, and the directionality in the photon emission is 
directly connected to the momentum conservation which, 
in turn, is a consequence of the constructive interference 
in the field emitted by each atom. Thus, the crucial is- 
sue in our scheme is to determine how this mapping is 
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modified in finite atomic ensembles under nonideal con- 
ditions. In the following we analyze such questions in 
detail, concentrating in the simplest case in which we 
have a single excitation with momentum |ko) in \sa) (i.e. 
our initial state is a W-like state) and thus we produce 
a single photon. We determine a function /(il), which is 
proportional to the probability density that the photon 
is emitted in the direction fl. In general, / = /coh + /inc; 
that is, it is the sum of a coherent contribution and an 
incoherent one. The later appears whenever the posi- 
tions of the particles fluctuate, /coh contains the forward 
scattering contribution, which is emitted in a cone with 
a width Af2 that decreases with the number of particles, 
/inc, on the contrary, describes isotropic light emission, 
thus, even when the light emitted in Ail is collected, the 
contribution /inc leads to a limitation in the efhciency of 
the setup. To quantify the error probability, we define 

^ /dri/incffi) , . 

As long as the number of excited atoms is small <C N, 
this analysis can be easily generalized to the emission of 
states with many photons (H]). One obtaines that the 
overall error is bounded by 1 — (1 — 

The emission pattern can be obtained by studying the 
Heisenberg equations of motion of the field operators. 
The calculation inolves the study of the decay of the 
atomic state under collective effects (see Appendix [C| . 
To simplify our analysis we ignore the dipole pattern, in 
which case we get: 

1 ^ 
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Tj are the coordinate operators of the atoms, and thus Eq. 
([S]) allows us to describe fluctuations in the position of the 
particles during the emission of light. In the following, 
we will show three different experimental set-ups where 
our scheme can be implemented. In order to analyze 
the performance in each of them, we first particularize 
the above formula to three different situations which are 
directly connected with those set-ups. We will focus on 
the angular width of the forward-scattering cone, AQ, 
which measures the collimation of the emitted photons, 
and the error probability, £, as figures of merit. Then 
we will introduce the possible implementations and will 
use those formulas to specify the conditions for them to 
correctly operate. 

(i) Fixed atomic positions. In the case of a square 
lattice of particles trapped in 3D (see Fig. [1] (b)), the 
emission pattern is given by 

f(n)^-Y\ sin^((fcLn^-(kL+ko)")dojV„/2) 

with Na the number of atoms in each direction. fo{^) 
has a series of diffraction peaks, which are reduced to a 



single one if lio < A/2. In this regime, the emission is cen- 
tered in a cone with no in the direction of Icl -I- kp. Note 
that for simultaneous energy and momentum conserva- 
tion condition |kL + ko| — ki^ has to be fulfilled. Since 
the positions of the atoms do not fluctuate, /o has only 
a coherent contribution {£ — 0), and the only limitation 
for the effiency of the setup is the width of the emission 
cone, which scales in 3D fike A93d ~ l/(A^^/^fcLrfo)- In 
the case of a chain of atoms (ID) momentum is conserved 
only along the direction of the chain. Photon emission 
can be still directed efficiently along the axis of the chain, 
in a cone whose width scales like A6'id ~ 1/^/Nki^do. 

(ii) Fluctuating atomic positions. Let us consider a 
lattice of atoms at temperature T, trapped by indepen- 
dent harmonic potentials. The emission pattern is now 
the sum the of two contributions, 

ho^^m = h{n)gTm, /inc(r!) = l-gT(f^). (7) 

gT{^) = e"((''Lnf!-kL)«T)"^ and is the vector whose 
components are the size of the position fluctuations in 
each spatial direction, (^f^)^ = x'^{l + 2ny), with Xq the 
size of the ground state in harmonic potential, and 
the number of motion quanta at T. Light scattered into 
/inc represents an important fraction whenever ^ do- 
In this case, the emission of light is centered around kL, 
since the uncertainty in the position of the particles av- 
erages out the intial linear momentum ko . The scaling of 
£ in this regime strongly depends on the dimensionality 
of the system. In particular, in the case of a chain of 
atoms, £iD = d^jX, whereas in the square 3D lattice, we 
get £3D « 12.6(do/A)27V-i/3. 

(Hi) Statistical distribution of particles. Consider an 
ensemble of atoms (see Fig. [T] (c)), which move inside 
a square box of size L, such that their motion is faster 
than their radiative decay, that is, their average velocity 
V is such that vL ^ F, with F the emission rate. This 
situation can be described by assuming that the atoms 
are in a statistical distribution with equal probability to 
be at any point in the box. The situation is thus similar 
to that of a thermal state, 

/coh(r!) =iVgbox(r!), /inc(^^) = l-gbox(f^), (8) 
where 

5box(f^)= n sinc^ (L(fcin^ - kg)) . (9) 

Oi—x.y,z 

Defining the average distance between particles like do = 
LN^'^/'^, we find the same scalings of AO as in case (i), 
and of £, as in case {ii). Trapping schemes for atomic 
ensembles are simpler to realize but face the difficulty 
that conditions for the directionality of photon emission 
are more stringent. In the case of a lattice of particles at 
fixed positions, forward-scattering is ensured whenever 
condition do < A/2 is fulfilled. On the contrary, in the 
case of atomic ensembles, the incoherent contribution /inc 
has to be small enough such that £ ^ 1, which implies 
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FIG. 1: (a) Level configuration for the release of atomic en- 
tangled states in photonic channels, (b) Release of a collec- 
tive state with linear momentum ko, that has been gener- 
ated in a lattice of atoms, (c) Emission of photons from an 
atomic ensemble, which consists of an incoherent contribu- 
tion (isotropic), and a coherent one in the forward-scattering 
direction. 



do ^ A in ID, or, alternatively, a number of particles 
large enough in 3D. 

Now we introduce three experimental set-ups where 
our scheme can be implemented. In the Appendix we 
show how to create the atomic states that we are consid- 
ering here. 

Trapped ions. This system is ideally suited to create 
collective states like ([1]), as was demonstrated recently 
in ref. fld\. Most usually ions are arranged in chains, 
such that we deal with the ID situation discussed above. 
Even though trapped ions are not equally spaced, under 
the condition dp < A/2, with do the average distance, 
we still get light emission in the forward-scattering cone 
only, see Fig. [2] Considering two different internal levels, 
which can correspond to different states in an hyperfine 
multiplet, states such as those defined by Eq. ^ can be 
created by a number of quantum operations that scales 
linearly with the number of ions A'^ (see lAip . For exam- 
ple, the state l/y2(|0,2kL) + |2ki,0)), would emit two 
photons in the forward and backward directions along the 
chain axis, entangled in polarization. The main difficulty 
for the implementation of this idea with ions lies on the 
fact that ion-ion distances are usually in the range of a 
few fim, and thus condition do < X/2 is not fulfilled when 
considering optical wavelengths. A way out of this prob- 
lem is to use optical transitions which lie in the range of 
A > 5 Mm, which can be found in ions such as Hg+, Ba+, 
or Yb+ dl. 

Cold atoms in optical lattices. By using optical lat- 
tices we fulfill the need of placing atoms at interparticle 
distances comparable to optical wavelengths, since po- 
tential wells in a standing-wave are indeed separated by 
do = Asw/2, with Asw, the wavelength of the counterprop- 
agating lasers. By using an optical transition such that 
A > Asw, we are in the regime in which light emission 
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FIG. 2: Probability of photon emission from an ion chain with 
= 30 ions initially in a W-state. The blue line corresponds 
to a chain with equally spaced ions with two diffraction peaks. 
Black and red lines corresponds to an ion Coulomb chain, in 
which ions are in an overall trapping potential and thus are 
not equally spaced. However, in the case that the average dis- 
tance, do is small enough, light is also preferentially emitted 
in the forward-scattering direction. 

is focused into a single Bragg peak. Although one could 
think of peforming quantum gates between ultracold neu- 
tral atoms to generate collective atomic states j22i], this 
procedure faces the difficulties of quantum computation 
in this system, like for example how to achieve single 
atom addressability. More efficiently, one could avoid 
the use of quantum gates by using the dipole-blockade 
mechanism with Rydberg atoms, which allows us to gen- 
erate W-states, as well as states which emit Fock states 
with a number M of photons ^IS'l (see Appendix [3]) . 

Atomic ensembles at room temperature. The very same 
techniques which can be applied to Rydberg atoms in an 
optical lattice can also be used in the case of hot en- 
sembles. On the one hand, this setup has the advantage 
that atoms do not need to be cooled and placed in an 
optical lattice. On the other hand, it can be described 
by a statistical distribution of particles, and thus suffers 
from the fact that high efficiency in the release of pho- 
tons is achieved under more severe conditions of particle 
density and atom number, as discussed above. However, 
densities which are high enough to fulfill the requirement 
£ ^ 1 have been recently reported in [23 |. 

In conclusion, we have proposed to use current tech- 
niques for quantum engineering to generate atomic multi- 
partite entangled states which can be efficiently mapped 
into photonic states. Our proposal relies on the release of 
spin-wave like excitations into a given spatial direction 
by means of interference effect, and can be implemented 
with trapped ions, atoms in optical lattices, and atomic 
ensembles at room temperature. 

This work was supported by E.U. projects (SCALA 
and CONQUEST), and the Deutsche Forschungsgemcin- 
schaft. 



4 



APPENDIX A: CREATION OF COLLECTIVE 
ATOMIC STATES IN A CHAIN OF ATOMS 

Entangled states of the form ([T]) and their hnear com- 
binations can be generated in a chain of particles, for 
example, of trapped ions, by means of a limited number 
of quantum operations. To demonstrate this, we first 
show that they can be written as Matrix Product States 
with a small bond dimension D, i.e. they can be written 
as 

\^)^ ^ ($F|y[-j...y[l]|$i)|ii)...M, (Al) 

ii,...,ijv 

In (jAip . the indices ij = g,Sa,Sb, and Vj^j are D x D 
matrices acting on an auxiliary D-dimensional Hilbert 
space. D is given by the number of states which appear 
in the singular value decomposition (s.v.d.) of j^*) at 
any site in the chain ^] . As it is shown in [2^ , the state 
(jAip can be prepared by performing N gates which act 
on [log2 £*] + 1 qubits. Thus, as long as D is independent 
of N, the number of gates to be applied scales linearly 
with the total number of atoms. 

To evaluate D, consider first the case of a state like ([T]) 
with atoms excited in level Sa only, and a partition of the 
chain in two parts L and R. We get tIq + 1 states in the 
s.v.d. with respect to this partition, which correspond to 
states with a number of excited atoms in part L, ranging 
from to This result is easily generalized to a linear 
combination of M states of the form ([T]), in which case we 
get D = M{na + l){ni, + 1). For example, an entangled 
state of the form: 

I*) = i= (|k"" = l, q"'=l) + |q"" = l, k"'=l)) , (A2) 
has D = 8. 



APPENDIX B: QUANTUM STATE 
ENGINEERING WITH RYDBERG BLOCKADE 

Interactions between excited atomic states, like those 
that take place in Rydberg atoms, can be used to the cre- 
ate the states defined by Eq. (|A2p . This can be achieved 
in a single experimental step, without the need for quan- 
tum gates, if the proper configuration of atomic inter- 
actions is chosen. As an example, consider the 3 level 
configuration shown in Fig. [T] (a), and interactions be- 
tween excited states such that atoms in levels |sq), \sb), 
interact strongly only if they are in the same excited 
state, that is, Uaa — Ubb — U, but Uab — 0. We ap- 
ply two lasers with wavectors ki^2 and Rabi frequencies 
r2i^2, detuned with respect to the \g) - \sa^b) transition, 
such that Al = — A2 = A. If condition Ai^2 ^ ^^1,2 
is fulfilled, then the lasers induce a two-photon transi- 
tion with Rabi frequency f2eff = riifi2/A. Furthermore, 
if r^off ^ t/, states with two atoms in the same excited 
state are not populated. Under these conditions there 



are two possible excitation channels, depicted in Fig. [3l 
which give rise to the linear combination (|A2p . 




FIG. 3: Lasers and level configuration for the creation of 
atomic entangled states which emit pairs of photons entan- 
gled in polarization. 



APPENDIX C: CALCULATION OF THE 
PHOTON DISTRIBUTION 

We consider for simplicity the lambda configuration 
depicted in Fig. [1] (a), considering a single excited state 
|s), and a single auxiliary level |e). The interaction of 
the quantized electromagnetic field with the ensemble of 
atoms, after the adiabatic elimination of level |e), is de- 
scribed by 

Hi = (^T]ak,Ae'('^-'^^>^+'"^* + h.c}j , 

^■-■^ = ^^^^ 

aj refers to the \g) - \s) atomic transition, fi^ and kL 
are the Rabi frequency and wave-vector of the classical 
field, respectively, ujk is the photon energy, ek,A are the 
polarization vectors, and dge is the dipole matrix element 
for the \g) - |e) transition. 

The probability of photon emission is proportional to 
the the diagonal elements of the one-photon density ma- 
trix, which are obtained from the Heisenberg equation of 
motion for the field operators, 

1 

^^^-.(k-k.)(r.-r,)^^^t(^^y^.(^^))_(C2) 
ij 

Since we are interested in the conditions for momentum 
conservation due to interference effects, we consider the 



5 



following atomic initial state, 

|ko)=<|0), < = ^Ee-'''-°at. (C3) 

The emission pattern depends thus on the two-time 
atomic correlation function, which in turns can be eval- 
uated by means of a master equation which describes 
the decay of the atomic levels. In the case of the ini- 
tial atomic state kg (jC3p . fixed atom positions, and ne- 
glecting boundary effects, this correlation function can 
be evaluated exactly, 

(aj(ri)a,(r2)) = e-^^'^^^+^^^/^e^'^°'^'^°-'^°\ (C4) 

where we have neglected an energy shift due to dipole- 
dipole interactions. Integrating (|C2[) over the absolute 
value of k yields the probability of photon emission, 

i{n) = i{n)f{n). (C5) 

/(O) is the dipole pattern, 

m) = ^^^{l-in,,nnf), (C6) 

where T is the single atom radiative decay rate, is 
the collective decay rate, ngg is the unit vector of the 
atomic transition, and nn is a unit vector in the direc- 
tion defined by the solid angle fl. The factor /(O) in 
describes the interference between the emission from dif- 
ferent atoms, and is given by Eq. ([5]). 

Below we deduce the master equation which leads to 
(|C4p and we sketch its solution in the case of collective 
states with a single excited atom. 

APPENDIX D: MASTER EQUATION 

The master equation for the reduced density matrix 
of the internal levels, which describes the ratiative decay 
of a set of atoms under the coupling to the quantized 



radiation field given by Eq. (|C1|I . is 

dtp = '^iP^l ~ '^l^iP ~ P'^i'^i) 

where the coupling constants depend on 

/>oo 

Jo 

in the following way: 

3(J,,) = ic,,. (D3) 
The master equation (jDl[) can be solved for the particular 
case of an initial state (jCSp by noticing that the evolu- 
tion of p is closed within the subspace spanned by the 
states |ko), |0). This fact can be easily proved by direct 
substitution of p(0) = |ko)(ko| in Eq. (jDl[) . which yields 
the following evolution for the atomic density matrix: 

p{t) = e-i^''o*|ko)(ko| + (1 - e-r^o*) |0)(0|, (D4) 

where the collective decay rate is just the Fourier 
transform of the coupling constants in the master equa- 
tion, 

ruo = Er..e*°(-"-"^ (D5) 

j 

A similar result holds for nondiagonal elements of p{t), 
which together with the quantum regression theorem 
yields the evolution of the atomic correlation function 
(El. 
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